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ABSTRACT 

We have evaluated the energetics of 38 solar eruptive events observed by a 
variety of spacecraft instruments between February 2002 and December 2006, as 
accurately as the observations allow. The measured energetic components in- 
clude: (1) the radiated energy in the GOES 1 - 8 A band; (2) the total energy 
radiated from the soft X-ray (SXR) emitting plasma; (3) the peak energy in the 
SXR-emitting plasma; (4) the bolometric radiated energy over the full duration of 
the event; (5) the energy in flare-accelerated electrons above 20 keV and in flare- 
accelerated ions above 1 MeV; (6) the kinetic and potential energies of the coronal 
mass ejection (CME); (7) the energy in solar energetic particles (SEPs) observed 
in interplanetary space; and (8) the amount of free (nonpotential) magnetic en- 
ergy estimated to be available in the pertinent active region. Major conclusions 
include: (1) the energy radiated by the SXR-emitting plasma exceeds, by about 
half an order of magnitude, the peak energy content of the thermal plasma that 
produces this radiation; (2) the energy content in flare-accelerated electrons and 
ions is sufficient to supply the bolometric energy radiated across all wavelengths 
throughout the event; (3) the energy contents of flare-accelerated electrons and 
ions are comparable; (4) the energy in SEPs is typically a few percent of the 
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CME kinetic energy (measured in the rest frame of the solar wind); and (5) the 
available magnetic energy is sufficient to power the CME, the flare-accelerated 
particles, and the hot thermal plasma. 

Subject headings: Sun: activity - Sun: coronal mass ejections - Sun: flares - 
Sun: particle emission - Sun: X-rays, gamma rays 



1. Introduction 



Solar eruptive events (SEEs), which are comprised of flares and associated coronal mass 
ejections (CMEs), are the most energetic occurrences in the solar system. Over a period 
of tens of seconds to minutes, they can convert upwards of 10 32 ergs of energy carried in 
non-potential, current-carrying magnetic fields into accelerated particles, heated plasma, and 
ejected solar material. 

While the overall energy involved in a large SEE is not in serious doubt, its partition 
amongst its component parts has so far been estimated only for a few events. In this 
paper, we provide the first statistical analysis of energy partition throughout the various 
manifestations of an SEE, for thirty-eight large events. We provide this information not 
only to establish "typical" ratios of the energy in various components of the event, but also 
to provide some idea of the range over which such ratios extend, and we especially point out 
events in which the strength of one component or another appears to lie outside the norm. 
We offer this analysis with the goal of providing useful constraints for modelers of the energy 
release process(es) involved. 

This paper grew out of the energetics working group at the meeting on "Solar Activity 
during the onset of Solar Cycle 24" held in Napa, CA, from December 8 - 12, 2008. It is a 
continuation of the wor k begun at the Tao s ACE /RHESSI/WIND joint workshop in 2003 
that led to the works of lEmslie et al.l ( 120041 . 120051 ) . These papers provided the first detailed 
analysis of most of the components of two well-observed SEEs (the GOES XI. 5 event on 
2002 April 21 and the X4.8 flare event of 2002 July 23), including the energies in thermal 
plasma, flare-accelerat e d elec trons and ions, associated CME, and solar energetic particles 
(SEPs). lEmslie et al.l ( 120041 ) showed that, for the two events in question, the energy in 



the magnetic field was sufficient to power the thermal soft X-ray (SXR) emitting plasma, 
the flare-accelerated electrons and protons, and the kinetic energy in the CME, and they 
also provided order-of-magn itude estimates of the partition of the energy amongst these 
components. Subsequently, lEmslie et al.l (120051 ) also considered the energy in the optical 
and EUV continua, and they cautioned that, due to the transfer of one energy component 
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to another (e.g., flare-accelerated electrons — > thermal plasma — > SXR emission), care must 
be taken in summing energetic components to arrive at the total energy rele ased in an SEE. 
The present paper is also motivated by the work of iMewaldt et al.l (l2008al ). which was the 
first to address the ratio of two energetic components (the CME energy in the rest frame of 
the solar wind and the energy in SEPs) for a statistically significant number of well-observed 
events. 

The basic objective of the paper is to conduct a statistical study of the energy partition 
into different components for many of the larger SEEs observed during the previous maxi- 
mum of solar activity, particularly during the period February 2002 through December 2006, 
the first fi ve years of obse rvations by the Ramaty High Energy Solar Spectroscopic Imager 
(RHESSI; iLin et al.ll2002[ ). The intent is to apply previously proven tec hniques to det e rmin- 
ing t he global energetics of many more events than the two studied by lEmslie et al.l ( 120041 
20051 ). and, where possible, to apply new techniques to improve the energy estimates. 



Our energy estimates come from a wide variety of observations: CME ki netic and poten 



tial e nergies from the Large Angle and Spectrometric COronagraph (LASCO; lBrueckner et al. 



19951 ) instrument on the Solar and Heliospheric Observatory (SoHO); energy in flare-accelerated 
charged particles inferred from the hard X-rays and gamma-rays observed by RHESSI; en- 
ergy contained in the SXR-emitting hot plasma from the Geostationary Operational Envi- 
ronmental Satellites (GOES) and RHESSI; energy in SEPs from the suite of instruments on 
the Advanced Composition Explorer (ACE) and from GOES, SoHO, the Solar Anomalous 
and Magnetospheric Particle EXplorer (SAMPEX), and the Solar TErrestrial RElations Ob- 
servatory (STEREO); a nd total radiated energy from the Total Irradiance Monitor (TIM, 
Kopp fc Lawrence] 120051 ) on the SOlar Radiation and Climate Experiment (SORCE). For 
weaker events, or where total irradiance measurements are not available, the Flare Irradi- 



ance Spectral Model (FISM; IChamberlin et al.ll2007l . 120081 ) was used to provide estimates of 
the bolometric output of a flare based on other measurements. 

In Section [2] we present the events studied and review the techniques used to estimate the 
different component energies of each event. In Section [3] we present a series of scatter plots 
of one energy component against another. While the uncertainties on the individual energy 
estimates are typically large (often an order of magnitude or greater), these scatter plots, 
because of the relatively large number of events they contain, nevertheless allow some general 
conclusions to be reached (Section H]) about how the energy is partitioned. The spread in 
the values for the different energy components also gives an idea of the uncertainties in the 
measured parameters and the range of flare intensities of the selected events. These plots 
also allow for the identification of a few "outlier" events (Section [5]) that indicate either 
larger measurement uncertainties or distinctly different energy partitioning for those events. 
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We summarize the results in Section [6], which also provides suggestions for future work. 



2. Component Energies of the Solar Eruptive Events 

The events studied are listed in Table [TJ They include the largest SEP events observed 
after February 2002, when RHESSI was launched, excluding those events beyond the West 
limb (for which no reliable active region identification can be made) and those events located 
from E60° to E90° (for which the evaluation of the SE P energy is h i ghly u ncertain - see 



Section 1277)) . They also include the two events studied by lEmslie et al.l ( 120041 ). which appear 



as Events #2 and 6 in Table [TJ Additional events include all flares for which RH ESSI 



detected significant (> 4cr) emission in the 2.223 MeV neutron-capture gamma-ray line (jShih 



20091 ; IShih et al.ll2009l ). This, plus the inclusion of an intriguing behind-the-limb event with 
a strong CME on 2002 July 20 (Event #5 in Tabled]), resulted in a total of 38 events for 
study. As permitted by the available data, estimates were made of the following energies for 
each of the 38 events: 

1. Radiated energy in the GOES 1 - 8 A band; 

2. Total radiated energy from the SXR-emitting plasma; 

3. Total (bolometric) radiated output; 

4. Peak thermal energy of the SXR-emitting plasma; 

5. Energy in flare-accelerated electrons; 

6. Energy in flare-accelerated ions; 

7. CME kinetic energy in the rest frame of the Sun; 

8. CME kinetic energy in the solar-wind rest frame; 

9. CME gravitational potential energy; 

10. Energy in SEPs; and 

11. Free (nonpotential) magnetic energy in the active region. 



It is important to keep in mind the differences among the first four items on this list. 
They are all related, but are included separately since they can each be estimated indepen- 
dently, and indeed relatively straightforwardly, from the available measurements, and since 
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collectively they provide significant information on the thermal energy of each flare, how 
it is distributed in temperature, and when it is generated and released. Further details on 
these four items are provided in subsections 12.11 to 12.31 Broadly speaking, the first item 
is the energy radiated in the narrow GOES band from 1 - 8 A, obtained directly from 
background-subtracted data (Section l2.ip . The second is the energy radiated over all wave- 
lengths (including the 1 - 8 A band) from the hot SXR-emitting plasma, and is a quantity 
inferred from the plasma parameters (emission measure and temperature) revealed by the 
GOES 1 - 8 A measurements. The third item is the total energy radiated over all wavelengths 
from all components of the flare at all temperatures (including that from the SXR-emitting 
plasma); in some cases this is directly observed and in some cases inferred from modeling of 
the emission in select wavelength ranges - see Section 12.21 The fourth item does not specify 
a radiated energy at all, but rather the peak thermal energy content of the hot SXR-emitting 
plasma; this quantity is inferred (Section |2.3|) from the pa rameters of spectroscopic fits to 



RHESSI data. It is important to realize (jEmslie et al.ll2005l ) that these four components are 



not separate flare energy components and therefore cannot be summed together to obtain a 
total flare thermal energy. 

We have not evaluated energy losses from the SXR-emitting plasma by thermal con- 
duction. However, we would note that conductive transfer of thermal energy from hot SXR- 
emitting plasma into the relatively cool chromospheric plasma will generally result in the 
thermal energy content of the hot SXR-emitting plasma being quite effectively transported 
to, and ultimately radiated away by, such relatively cool plasma, one contribution to the 
total (bolometric) radiated output - the third item on the list. Further, for the events con- 
sidered here there is little observational data available on the energy contained in turbulence 
and directed mass motions of thermal plasma, components that may well c ontain energies 



comp arable to the thermal energy of the SXR-emitting plasma (see, e.g., iDoschek et al. 



1992|). 



The data on all the component energies are summarized in Table [TJ GOES SXR data 
are available for all the events. However, because of missing or inadequate data or limited 
instrument sensitivities, definitive energy estimates for all the energy components listed 
above are available for only six events (Events #13, 14, 20, 23, 25, and 38). As mentioned 
above, Event #5 was located behind the limb; thus only a "plausible" magnetic energy 
estimate (not included in Table [I]) could be obtained from observations of the most likely 
responsible active region once it had moved onto the solar disk, and the listed radiated 
energies are lower limits. 
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Table 1. Event List with Component Energies (x 10 30 ergs) 



No. 


Date* 


Time** 


Class 


SXR 1 


T-rad 2 


Bol 3 


Peak 4 


Elec 5 


Ion 6 


KE 7 


SW 8 


PE 9 


SEP 10 


Mag 


1 


02/02/20 


05:52 


M5.1 


0.043 


1.2 


13 








17 


5.6 


6.3 


0.13 


1200 


2 


02/04/21 


00:43 


XI. 5 


1.2 


38 


150 


13 


20 




230 


160 


5.0 


23 


660 


3 


02/05/22 


03:18 


C5.0 


0.048 


5.6 


9 








84 


45 


10 


2.7 


260 


1 


02/07/15 


19:59 


X3.0 


0.31 


6.4 


44 


>2.2 


>3.6 




160 


76 


10 


3.8 


1500 


5 


02/07/20 


21:04 


X3.3t 


>1.5 


>26 


>210 








260 


170 








6 


02/07/23 


00:18 


X4.8 


1.2 


19 


150 


2.5 


32 


39 


260 


150 


20 


<30 


2000 


7 


02/08/24 


00:49 


X3.1 


1.1 


21 


160 


5.9 


11 




210 


130 


16 


3.9 


2500 


8 


02/11/09 


13:08 


M4.6 


0.11 


5.0 


8 


1.3 


60 




180 


110 


20 


0.51 


550 


9 


03/05/27 


22:56 


XI. 4 


0.16 


3.6 


16 


2.8 


7.4 


0.19 










260 


10 


03/06/17 


22:27 


M6.9 


0.21 


4.6 


17 


2.4 


1.6 


6.7 










110 


11 


03/10/26 


17:21 


XI. 2 


1.2 


31 


88 








240 


1.30 


32 


0.75 


1700 


12 


03/10/28 


09:51 


X17 


4.4 


68 


362* 


>19 


>56 


>190 


1200 


850 


6.3 


43 


2900 


13 


03/10/29 


20:37 


X10 


1.9 


31 


137* 


11 


110 


30 


340 


220 


25 


9.7 


2900 


11 


03/11/02 


17:03 


X8.3 


1.8 


24 


1.30 


9.3 


1.30 


68 


270 


200 


10 


9.3 


2800 


15 


03/11/03 


09:43 


X3.9 


1.1 


17 


97 


2.4 


120 


3.1 










780 


16 


03/11/04 


19:29 


X28 


4.8 


72 


426* 


>3.1 


>21 




610 


410 


25 


5.3 


2800 


17 


04/07/15 


18:15 


XI. 6 


0.16 


1.1 


8 


0.93 


42 


<0.1 










820 


18 


04/07/25 


05:39 


M7.1 


0.069 


1.3 


10 














2.9 


2300 


19 


04/11/07 


15:42 


X2.0 


0.32 


5.0 


56 


.3.0 


43 




220 


130 


25 


4.2 


610 


20 


04/11/10 


01:59 


X2.5 


0.32 


7.7 


15 


2.0 


20 


3.4 


230 


180 


16 


2.4 


610 


21 


05/01/15 


00:22 


XI. 2 


0.23 


4.7 


23 


5.0 


.32 












1500 


22 


05/01/15 


22:25 


X2.6 


1.3 


22 


78 


7.1 


63 


15 


730 


540 






1600 


23 


05/01/17 


06:59 


X3.8 


1.8 


34 


150 


17 


48 


13 


1000 


730 


50 


11 


1600 


24 


05/01/19 


08:03 


XI. 3 


0.43 


7.0 


54 


5.9 


82 


29 










1600 


25 


05/01/20 


06:36 


X7.1 


2.9 


43 


150 


10 


25 


120 


15 - 79 


7.8 - 61 


2.0 


7.8 


1600 


26 


05/05/13 


16:13 


M8.0 


0.44 


14 


49 


3.1 


13 




39 


22 


4.0 


7.3 


400 


27 


05/07/14 


10:16 


XI. 2 


0.64 


12 


87 


4.3 


24 




100 


66 


6.3 


2.9 


310 


28 


05/07/27 


04:33 


M3.7 


0.16 


4.5 


.30 


1.3 


12 




100 


62 


10 




310 


29 


05/08/22 


16:46 


M5.6 


0.34 


9.8 


35 


.3.2 


6.3 




110 


76 


10 


6.4 


390 


30 


05/08/25 


04:31 


M6.4 


0.050 


1.2 


11 


1.1 


16 


<1.9 










110 


31 


05/09/07 


17:17 


X17 


4.9 


68 


322* 


>5.6 


>10 


>0.7 










1400 


32 


05/09/09 


19:13 


X6.2 


3.1 


41 


250 


>7.9 


>120 


>1.7 










1300 


33 


05/09/10 


21:30 


X2.1 


0.99 


17 


82 


6.0 


13 


1.0 










1300 


34 


05/09/13 


19:19 


XI. 5 


1.1 


25 


85 








330 


200 


32 


3.0 


1400 


35 


05/09/13 


23:15 


XI. 7 


0.23 


4.7 


21 


2.3 


32 


<0.1 










1400 


36 


06/12/05 


10:18 


X9.0 


1.4 


19 


92 


>5.1 


>360 


>4.5 










400 


37 


06/12/06 


18:29 


X6.5 


1.1 


18 


59* 


6.8 


40 


36 










110 


.38 


06/12/13 


02:14 


X3.0 


1.1 


17 


75 


4.8 


1.3 


14 


74 


41 


6.3 


.3.2 


570 



*In yy/mm/dd format. ** GOES start time (UT). 

Radiated energy in the GOES 1 - 8 A band. 2 Total radiated energy from the SXR-emitting plasma. 
3 Bolometric radiated energy. 4 Peak thermal energy of the SXR-emitting plasma. 
5 Energy in flare-accelerated electrons. e Energy in flare-accelerated ions. 

7 CME kinetic energy in the rest frame of the Sun. 8 CME kinetic energy in solar-wind rest frame. 

9 CME gravitational potential energy. 10 Energy in SEPs. 11 Nonpotential magnetic energy in the active region. 

*Behind-the-limb event. fBolometric irradiance directly measured with TIM - see Table [2] 
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2.1. Radiated Energy from Hot Plasma 



For each event, we estimated the time-integrated SXR and total radiated energies from 
the hot SXR-emitting plasma - the columns labeled 'SXR' and 'T-rad', respectively, in 
Table [TJ Fluxes in W m~ 2 for the 1-8 A and 0.5-4 A bands are provided by one of NOAA 
GOES satellites every 3 s. The total emission in the 1-8 A band ('SXR') is obtained simply 
by summing the background-subtracted fluxes over the duration of the flare, from the GOES 
start time (given by NOAA and listed in Table [1]) to the time when the flux had decreased 
to 10% of the peak value. The background that was subtracted was taken as the lowest flux 
in the hour or so before and/or after the flare. 

To calculate the radiated energy from the hot plasma, we used the measured GOES SXR 
fluxes in a manner similar to that described in lEmslie et al.l ( 12004] ). specifically using the IDL 
GOES Workbench available in SolarSoftware (SSW). This allows us to obtain a consistent 
set of values for all events since GOES, unlike RHESSI, has full coverage for all events. This 
calculation assumes that the hot plasma at any given time is isothermal; the temperature 
and emi ssion measure are calculated from the two-channel GOES data using the relations 
given by lWhite et al.l (120051 ) . Using the thus- inferred emission measure and temperature, and 
the optically th i n rad iation loss rate vs. temperature function (for coronal abun dances and 



1997 



Mazzotta et al.l (119981 ) ionization equilibria) taken from the CHIANTI database (jPere et al. 



2009 



j), we used the IDL procedure radJoss, available in SSW, to calculate, for each 



3-s time interval, the energy radiated from the SXR-emitting plasma over all wavelengths. 
Finally, we summed the radiated energies over the duration of the flare (from the GOES start 
time until the 1-8 A flux decreased to 10% of its peak value) to produce the total radiated 
energy given in the column labeled 'T-rad' in Table HJ Significant energy could be radiated 
af ter this nominal end of t he flare, particu larly if there is a "second phase" that, according 



to 



Woods et al.l (120111 ) and ISu et al.l (120111 ) , can release an amount of energy that is similar 



to that released in the initial phase. Generally, however, the values quoted in Table [TJ should 
include more than 50% of the energy radiated by the SXR-emitting plasma. 



2.2. Bolometric Irradiance 



Estimates of the bolometric irradiance, the total energy radiated from the flare inte- 
grated across the entire solar spectrum, for each of the events are provided in the column 
labeled 'Bol' in Table [TJ For five of the events listed in Tab le [TJ the bolometric irrad iance was 
measured directly by the Total Irradiance Monitor (TIM, iKopp fc Lawrence! 120051 ) onboard 
SORCE as an increase in the total solar irradiance above the (highly variable) pre- and 
post flare background levels. Total flare irradiance values were reported in this manner for 
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Events #12, 13, 16, and 31 by iWoods et al.l (120061 ) . and the bolometric irradiance for event 
#37 will be reported by Moore et al. (2012, in preparation). Both previously published 
values and revised estimates made for this paper are listed in Table [2J A final correction 
factor (see Table was then applied to allow for limb-darkening absorption when the path 
to the observer becomes optically thick at some wave lengths; the value of this factor can be 
up to ~ 3.0 - see Equation (2) in I Woods et al.l (120061 ). 



To complement these direct measurements and so provide a consistent set of bolometric 
values f or all of the events in Table [I I estimates from the Flare Irradiance Spectral Model 
(FISM; Chamberlin et al. 2007 , 20081 ) were used, with various assumptions and corrections 
as described below. FISM is an empirical model that provides estimates of the total amount 
of solar radiated energy over a broad wavelength range from 1-1900 A and over a wide 
range of time scales from seconds to y ears. It uses measu rements in this wavelength range 
from the Solar EUV Experiment (SEE; IvVoods et aUbood ) on the Thermosphere Ionosphere 
Mesosphere Energetics and Dynamics (TIMED) satellite an d the SOLar-STellar Irradiance 
Comparison Experiment (SOLSTICE; iRottman et al.lll993l ) on the Upper Atmosphere Re- 
search Satellite (UARS). 

For the relatively rapid GOES 1 - 8 A SXR flux variations that occur during a solar flare, 
different empirical factors appropriate to the rise and decay phases of the flare, respectively, 
are used to relate the SXR flux to the total radiated energy during those phases. Various 
daily proxies are also used to represent the more gradual variations in solar irradiance due 
to active region evolution, solar rotation, and the solar cycle. The daily pre-flare irradiance 
spectrum is subtracted from each value to get the radiated energy from the flare alone, and 
this is then integrated over the duration of the GOES flare to get the total radiated energy 
in erg cm -2 at the detector. Then, assuming uniform radiation over 2ir steradians, the total 
radiated energy from the flare in the 1-1900 A wavelength range can be calculated, with 
1-minute cadence. 

The 1-1900 A solar irradiance is converted to total radiated energy over all wavelengths 
(the bolometric irradiance) by multiplying by an empirical conversion factor of 2.42 ± 0.31, 
determined by comparing the 1-1900 A solar irradiance with the absolute bolometric intensity 
for the five flares for which the latter could be measured directly (see Table [2]). 

The uncertainties on the calculated values of the bolometric irradiance listed in Ta- 
bles [1] and [2] are made up of several parts. The most dominant uncertainty comes from the 
TIM measurements themselves, and is due to the variations in the total solar irradiance of 
the non-flaring Sun. Other contributions to the overall uncertainty are the errors on the 
FISM estimates of the 1-1900 A flux, the conversion from UV irradiance to total solar irra- 
diance, and the limb-darkening correction. The overall uncertainty on the calculated values 
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is ± 



^70% for those events that are near disk center and ± ~90% for the near limb events0 



2.3. Peak Thermal Energy Content of the Hot Plasma 



The peak thermal energy content of the SXR-emitting plasma (the column labeled 
'Peak' in Tabled]) was determined from RHESSI imaging spectroscopy data. We first fit the 
observed RHESSI hard X-ray spectra with the sum of a single-temperature Maxwellian plus 
the form expected from a double-power-law electron spectrum (Equation (J2J) in Section |2~4]) . 
The fit parameters appropriate to both thermal and nonthermal components were determined 
for each time interval using the forward-fitting method implemented in the OSPEX software 
package available in SSW. The temperatures and emission measures obtained from RHESSI 
in this way tend to agree closely with the corresponding values obtained from the standard 
GOES data analysis discussed in Section 12.11 However, somewhat higher temperatures can 
be obtained because of the RHESSI coverage to higher energies; indeed, superhot compo- 
nents with reported temperatures as high as ~50 MK may exist in some flare s and are not 
accur ately reflected in the GOES thermal analysis (ILin et al.lll98ll ; ICaspill2010t ICaspi &: Lin 
20101 ). The thermal function included the line plu s continuum componen ts determined using 



CHIANTI, again with coronal abundances and iMazzotta et al.l (Il998l ) ionization equilib- 
ria. From these fits, the average temperature T (K) and emission measure EM = J n 2 e dV 
(cm -3 ) of the thermal plasma were determined every 20 s throughout the flare. (Here n e is 
the electron density (cm -3 ) and V is the emitting volume (cm 3 ).) 

The thermal energy content Z7th of the plasma can be calculated from the expression 
(e.g, Ide Jager et al.lll9861 ) 



U th = 3n e kT fV ap ~ 4.14xlO- 16 T ov /£Mx/K P erg, 



(1) 



where k is Boltzmann's constant, / is the volumetric filling factor and V ap is the apparent 
volume of the SXR sourc e. Starting in 2003 , SXR images are also available from the GOES 
Soft X-ray Imagers (SXIs; lLemen et al.ll2004j; iPizzo et al.ll2005f) ; such images could be used to 



provide estimates of V ap , as described in lEmslie et al. l2004h . However, both for consistency 



with the earlier analysis of lEmslie et al.l ( 120041 ) (which analyzed events that occurred in 
2002, prior to the SXI deployment), and since the parameters EM and To in Equation (CQ) 
are deduced from RHESSI data, we have chosen to use estimates of K P that are deduced 



1 Because the conversion of the FISM radiated energy to bolometric energy is based on the five events 
measured directly with TIM, the bolometric energies for these five events derived from the FISM estimate 
differ (after correcting for limb darkening) by less than 40% from the directly measured values. 
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from RHESSI images, made using the 3<7 -clean metho d of iDennis &: Pernakl (120091 ) . We 
further take / to be unity , consis tent with Emslie et al.l (12004 ) and further justified by the 
recent work of iGuo et al.l (120121 ). who used hard X-ray imaging spectroscopy data of 22 
extended-loop events to derive a (logarithmic) mean filling factor / = 0.20 x/-=- 3.9 (la 
standard error). 

The peak energy values listed in Table [1] are the highest values of {7 t h obtained from 
this analysis, usually at or near the time of the peak GOES flux. 



2.4. Flare-Accelerated Electrons 

The energies in flare-accelerated electrons are listed in the column labeled 'Elec' in 
Table [TJ They were determined by using the OSPEX algorithm to fit a combined isothermal- 
plus-nonthermal function to the measured RHESSI spatially integrated X-ray spectra. The 
nonthermal component was assumed to be bremsstrahlung from energetic electrons with an 
injected spectrum Fq(Eq) (electrons cm" 2 s" 1 keV" 1 ) in the form of a broken power-law: 



Fq(E 



A < 





(E /E p ) 
(E /E p ) 




-<5i 
-S 2 



[E h jE l 



\<5 2 -<5i 



, Eq < -^min 

< E < E b 
, Eb < Eq < E max 
, E ma , x < Eq 



(2) 



The seven parameters of this model spectrum are the normalization parameter, A, the low 
and high energy cutoffs, E m i n and -E max , the break energy Eb, and the power-law indices Si 
and 82 below and above the break energy, respectively. The (arbitrary) value of the pivot 
energy E p was fixed at 50 keV. Also, the high energy cutoff E max was fixed at 30 MeV, an 
energy so high above the energy range of interest ( < 500 keV) that it has a negligible effect 
on the calculated X-ray spectrum and so is equivalent to having no high-energy cutoff at all. 

The OSPEX analysis uses a forward-fitting procedure that starts by dividing the flare 
into multiple time intervals - here we used 20 s intervals. For each interval, the function 
thick2 in SSW is used with a set of starting parameters for the electron spectrum (j2J) used 
to calculate the X-ray photon spectrum, assuming electron-ion bremsstrahlung in a thick 
target that is "cold" in the sense that the ambient electrons have a mean energy kT that 
is significantly lower than the lowest energy of the accelerated electrons. In general, con- 
sideration must also be given to the ionization state of the target, since the bremsstrahlung 



Brown 


1973; 


Kontar et al. 


2003) 



lower energy electrons that stop higher in the corona, we used parameters appropriate for a 
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fully-ionized plasma to calculate the total nonthermal energy . A mo re refined calcu lation is 



possible using the procedure outlined by iKontar et al.l (120021 ) and ISu et al.l ( 1201 ll ). but no 



significant difference is expected in the resulting total energy in electrons above E min . 

The resulting photon spectrum is then folded through the detector response matrix to 
generate a count-rate spectrum, which is added to the count-rate spectrum calculated for 
the thermal spectrum discussed in Section 12.11 Then, through an iterative procedure, we 
find best-fit values of the parameters describing the electron spectrum ([2]), by minimizing the 
X 2 statistic between the calculated and the measured background-subtracted count spectra. 
The total energy U e in electrons for a given event is then computed by integrating the best-fit 
electron energy spectrum above E m i n for each time interval and summing the results over 
the duration of the flare, resulting in the values listed in the column labeled 'Elec' in Tabled! 

In order to obtain the most reliable spectral fits to the RHESSI data and thus better 
evaluate the uncertainties in the calculated values of U e , we chose to use data from just one 
of RHESSFs nine detectors - Detector #4. This particular detector has good energy reso- 
lution and sensitivity, which allowed us to apply the most up-to-date corrections for energy 
resolution and calibration, photospheric albedo, pulse pile-up, and background subtraction 
that are available with the current analysis software. For the large events studied, the count 
rates were sufficiently high that selecting just a single detector did not seriously degrade 
the spectro scopic capability up to the photon energies required to determine the parameters 



of interest. iMilligan fc Dennis! ( 120091 ) have shown that similar best-fit parameter values are 
determined using different individual detectors (detectors 1, 3, 4, 5, 6, and 9 in their case), 
which leads to an estimate of the systematic uncertainties in the calculated total energy in 
electrons of ~20%. This is negligible compared to the uncertainty arising from the difficulty 
in establishing the value of the low energy cutoff energy E min , an uncertainty that arises 
because the thermal emission generally dominates the low-energy part of the X-ray photon 
spectrum up to energies where the effects of a cutoff in the electron spectrum might be de- 
tectable. We used the largest value of E m[n that still gave an acceptable fit (reduced x 2 — !)• 
As a result, the values of U e listed in Table [1] are lower limits to the energy in the nonther- 
mal electrons. Furthermore, because of the steep form of the electron spectra (5i <^4), these 
values are particularly sensitive to E min , so that the energies in flare-accelerated electrons 
could be up to an order of magnitude higher than those reported in Table [TJ 



2.5. Flare- Accelerated Ions 



The energies in flare-accelerated ions with energies above 1 MeV are listed in the column 
labeled 'Ion' in Table [TJ In order to provide a consistent set of values for as many events 
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as possible, the energies were estimated solely from RHESSI measurements of the fluence 
(time-integral of the flux, photons cm -2 ) in the 2.223 Me V neutron-c aptu re gamma-ray line . 
Our sample of events is primarily based on the studies of IShihl ( 120091 ) and lShih et al.l ( 120091 ). 
who analyzed RHESSI flares from 2002 to 2005 that had either 2.223 MeV line emission 
and/or >0.3 MeV electron bremsstrahlung continuum emission. Of those flares, energies 
are included only for those that have >2cr detections of that line, and further only as 4<r 
upper limits if below a 4cr detection. We also include three additional flares that occurred 
in 2006 (Events #36, 37 and 38). We chose a lower energy threshold of 1 MeV because 
the production of detected nuclear gamma-ray lines from elements such as 20 Ne begins at 
energies as low as ~3 MeV, and it is therefo re evident that the ion spectrum extends down 
to ~1 MeV, at least in a few large events (IRamaty et al.l Il995l ; iRamaty fc Mandzhavidze 



2000). The spectral shape is essentially unknown below 1 MeV. 



In order to estimate the energy in ions from the 2.223 MeV line fluences, the following 
steps were taken. The measured fluences of the line were first corrected f or attenuation in the 
solar atmosphere assuming a given depth of production of the photons (IHua fc Lingenfelter 
19871 ) and allowing for the flare position on the solar disk. The corrected fluence values 
were then con v erted to t he proton en ergy above 30 MeV using conversion factors given by 
Murphy et al.l (120071 ) and IShihl (120091 ) . The 30 MeV threshold was used at this stage in the 
analysis because the 2.223 MeV line is produced by ions with energies ^ 20 MeV nucleon -1 , 
so that the conversion factors are less dependent on the assumed power-law index of the 
proton spectrum. 

In order to estimate the energy in protons above 1 MeV, an extrapolation is required over 
one-and-a-half orders of magnitude in proton energy, so that the inferred energy above 1 MeV 
depends critically on the spectral index used in this extrapolation. For the largest RHESSI 
flares, where multiple types of ion-associated gamma-ray emission can be detected and fit 
simultaneously, the ion power- law spectral indices are found to be typically in the range 3-5, 
a range of indices consistent with that found in a study of flares obse rved by the Gamma-Ray 
Spectrometer on the Solar Maximum Mission ( Ramaty et ali 1996 ). Consequently, for the 
purposes of estimating the total energy in protons, we have assumed a power-law proton 
spectrum with a single spectral index of 4 that extends down to a lower cutoff energy of 
1 MeV. Because of the long "lever arm" associated with this extrapolation, an uncertainty 
in the spectral index of ±1 corresponds to an uncertainty in the total energy content above 
1 MeV of about ±1.5 orders of magnitude. 

Even under the assumption that the spectra for the various types of ions have the same 
spectral index and low-energy cutoff, the energy content will also depend on the accelerated 
particle composition. The ratio of the energy content in all ions (protons plus a-particles 
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and heavier nuclei) to the energy in protons can vary between ~2 and ~6; here we assume 
that the energy in flare-accelerated ions is three times the energy content in flare-accelerated 
protons. 

For a number of the events, the total energy content of ions listed in Table [1] is a lower 
limit because RHESSI did not see the complete time history as the result of spacecraft night 
or passage through the South Atlantic Anomaly (SAA). In addition, there can be other 
complications that affect the observation or interpretation of the neutron-capture line. The 
affected events are as follows: 



Event #12: RHESSI missed a significant fraction of the neutron-capture l ine emission. 



inclu ding the peak, as shown by observations of this flare by INTEGRAL fjKiener et al. 



2006|); 



Event #31: RHESSI observed only ~2 minutes of a significantly longer impulsive 
phase. Furthermore, the level of atmospheric attenuation is very uncertain due to this 
flare's large heliocentric angle, so we use a conservative angle of 80° to determine the 
correction factor; 

Events #32 and 37: RHESSI missed the peak of the impulsive emission, and thus 
possibly a significant fraction of the total emission; 

Event #36: RHESSI missed some fraction of the 2.223 MeV emission as it was just 
coming out of Earth shadow. RHESSI observations started at 10:31 UT and the GOES 
X-ray flare started at ~10:19 UT; 

Events #14, 15, 22, 33, and 38: RHESSI likely missed a small fraction of the neutron- 
capture line emission late in the flare. However, this missing energy is smaller than 
the other uncertainties in the energy estimates discussed above; 

Events #36, 37, and 38: By December 2006, RHESSFs detectors had reduced gamma- 
ray sensitivity resulting from accumulated radiation damage, a reduction that is diffi- 
cult to estimate. 



2.6. Coronal Mass Ejection 

The CME kinetic energies, both in the rest frame of the Sun and in that of the solar wind 
(for comparison with the SEP energies - see Section [3TB1) . are listed in the columns labeled 
'KE' and 'SW,' respectively, of Table [TJ The gravitational potential energies of the CMEs 
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are listed in the column labeled 'PE.' These C ME energies wer e estimated from calibrated 
LASCO images using the procedure detailed in IVourlidas et all (120101 l201ll ). 



Briefly, this procedure consists of the following steps. First, we selected two LASCO 
images, one containing the CME and the other taken before the event as close in time as 
possible to the flare with no disturbances or ejecta over the path of the subsequent CME. 
Next, the images were calibrated (in units of mean solar brightness) and the pre-event image 
was subtracted from the CME image. The excess brightness revealed by this subtracted 
image is due to Thompson scattering of photospheric radiation from the excess mass in the 
CME. This excess brightness can therefore be converted to excess mass of the CME under 
the usual assumptions that (1) all of the CME mass is co ncentrated on the plane of the sky, 



and (2) the CME material consists of 90% H and 10% He (IPoland et al.lll98ll : IVourlidas et al. 



2 OOOl . 120100 . We used the first assumption because the true three-dimensional distribution 
of the CME mass along the line of sight is unknown. The second assumption represents an 
"average" coronal composition, since we do not know the height at which the bulk of the 
CME material originates (other than that it is coronal). 

These assumptions together result in a lower limit for the mass. The uncertainty in 
the CME mass becomes more significant as the central angle and/or spread of a given CME 
departs significantly from the plane of the sky. T he mass underestimati on is about a factor 
of 2 for CMEs that are <40° from the sky plane (IVourlidas et al.ll201ol ). 



Other uncertainties in this procedure include exposure time variations between event 
and pre-event images, improper vignetting correction, solar rotation effects, and the presence 
of stars in the field of view. Fortunately, such uncertainties can be minimized to a level that 
is well below that of other factors through proper calibration and careful choice of event and 
pre-event images. 

After obtaining a series of excess masses of the CME as a function of time, we can 
compute both the total mass of the CME and the position and projected velocity, both for 
the leading edge and for the center of mass of the CME. From the mass m, position r, and 
velocity V we can straightforwardly estimate the total kinetic (Uk = \mV 2 ) and potential 
(£/$ = GMqUiIRq 1 — r -1 ]) energies (here G is the Newtonian gravitational constant and 
Mq,R q are the solar mass and radius, respectively). Thes e values are aga i n lowe r bounds 
since both the mass and the speed are projected quantities. IVourlidas et all (120101 ) estimate 
that, for CMEs that are far away from the sky plane and that have relatively small widths, 
the kinetic energy could be as much as a factor of eight times larger than the values derived 
above; similarly the potential energy could be as much as twice as large for such events. 
However, for the majority of events, the uncertainties on the quoted energies are within a 
factor of two. To obtain the kinetic energy in the solar wind rest frame (as an estimate of 
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the energy available for shock acceleration of SEPs; see Section [376]) . we simply subtracted 
400 km s -1 from the measured CME speed and recomputed the kinetic energy using the 
speed in this new reference frame. 



2.7. Solar Energetic Particles (SEPs) 



For the majority of the events studied, it is likely that the interplanetary SEPs in 
the events studied are accelerated by CME-driven shocks. (A possible exception is the 
200 2 February 20 even t, where particles directly accelerated in the flare could dominate; 



sec 



Chollet et al.ll2010l ). The energy content of the accelerated SEPs, particularly when 
compared to the kinetic energy of the CME in the solar wind rest frame, is therefore an 
important measure of the efficiency of SEP production by the CME. 

The energy content of SEPs that escape into interplanetary space has been estimated 
by measuring the energy spectra of electrons from ~0.035 to ~8 MeV, protons from ~0.05 
to ~400 MeV nucleon -1 , and abundant heavier ions from ~0.05 to ~100 MeV nucleon -1 . 
Estimates were made in a number of large events from Solar Cycle 23 using a combination of 
nine separate instruments. The proton spectra are based on data from the Ultra-Low Energy 
Isotope Spectrom eter (ULEIS ; 



Mason et al.l Il998l ). and the Electron, Proton, and Alpha 



Coo 


c et al. 


1993) 


1996 


) on NOAA' 



measured by the Solar Isotope Spectrometer (SIS; IStone et al.lll998l ) on ACE and by ULEIS. 



Also used for the 2006 even ts were two STEREO instruments , the Low Energy Telescop e 



fLET; lMewaldt et al.ll2008bf ) and High Energy Telescope (HET: Ivon Rosenvinee et al.ll2008r ). 



Electron measurements were provide d by ACE /EPAM, SAMPE X/PET, and by the Electron 



Proton Helium INstrument (EPHIN; iMuller-MeUin et al.lll995l ) instrument on SoHO. 



For eleven of these events, the energy spectra of H, He, and abundant heavi er ions were 



all fit with co mmon spectral forms that include the double-power-law function of iBand et al. 
( 119931 ) and the Ellison fc Ramatyl ( 119851 ) spectrum - a power- law with an exponential cutoff - 



Exam ples of energy sp e ctra a nd both functional forms are given in iMewaldt et al.l ( 12005 



2012 


) and 


Cohen et al. 


(2005) 



were fit and the contributions of He and heavier ions were estimated using element abun- 
dances measured for these events by ULEIS and SIS. The electron contribution was measured 
in each of the individual events using either EPAM and PET, or EPAM and EPHIN. 



For all of the fluence measurements described above, the instruments were located near- 
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Earth. As in lEmslie et al.l ( 12004] ). we used the measured near-Earth fluence spectra, typically 
integrated over 3 to 5 days, to estimat e the energy c m" 2 t hat escaped beyond 1 AU in the 



form of SEPs. To obtain this estimate, lEmslie et al.l (120041 ) corrected for the fact that SEPs 
can scatter back and forth across 1 AU (providing multiple opportunities to be measured) 
using correction factors based on simulations by Giacalone (personal c ommunication, 2002) . 
A similar approach was followed in analy zing the "Hallo w een" events (IMewaldt et al.l 120051 ) 
and in a subsequent survey of 17 events (lMewaldtll2006l ). IMewaldt et al.l ( ]2008a[ ) improved 
on these estimates in a study of 23 SEP events from 1997-2005 by correcting for the fact that 
SEPs also lose energy as they scatter on the diverging interplanetary magnetic field (IMF). 

For this w ork, we correct e d for both multiple 1-AU crossings and energy loss using new 
simulations by Chollet et al.l (2010 ) for fo ur species (H, He, O, and Fe) with a range of 
charge-to- mass ratios. Ichollet et al. ( 2010 ) considered scattering mean free paths A ranging 
from 0.01 to 1 AU, and also varied the radial and rigidity dependence of A. Surprisingly, the 
source energy required to account for the accelerated particles in these different scattering 
descriptions varied by less than a factor of ~2. This is apparently because the scattering and 
energy-loss processes compensate for each other - the more particles scatter the more often 
they cross 1 AU, but they also lose more energy in the process. In t his paper we have u sed 
their form of A derived from quasi-linear theory (see Equation (3) in I Chollet et al.ll2010l ). 



To relate the measured near-Earth values of MeV cm -2 to the integrated contribution 
of SEPs escaping through a 1-AU sphere surrounding the Sun, we nee d to know how SEPs 
from a given source location are distributed in longitude and latitude. lEmslie et al.l (120041 ) 
assumed that the SEP fluence at Earth falls off exponentially with e-folding separations 
of 35° for latitude, 45° for longitude in western events, and 25° for longitude in eastern 
events. Since then, lLario et al.l (120061 ) have measured the longitudinal distribution of SEPs 
using two- and three-spacecraft data from the two Helios spacecraft and the Interplanetary 
Monitoring Platform-8 (IMP-8). They adopted a Gaussian spatial distribution given by 
F = F a exp[— ($ — $ ) 2 /2<7 2 ], where $ is the longitude of the observer, $ is located 25°. 8 
east of the point of best solar wind connection for a 450 km s _1 solar wind (~W52°), and 
a = 38°. We use their result for the fluence of 4 - 13 MeV protons and we assume it also 
applies to latitude differences. By using this relation with the measured flare location and 
the near-Earth value for the escaping MeV cm -2 , we obtained the source energy required to 
supply SEPs escaping over a 1-AU sphere centered on the Sun. The results are tabulated 
in the 'SEP' column of Table [TJ Note that we have limited this study to SEP events with 
source regions ranging from E60 to W90 in longi tude. Beyond this range the Gaussian 
longitude distribution adopted by lLario et al.l ( 120061 ) drops off very rapidly and the longitude 
corrections become considerably greater and more uncertain. The typical uncertainty in SEP 
energy is conservatively estimated to be a factor of three. 
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2.8. Nonpotential Energy in the Active Region Magnetic Field 



It is commonly believed that the fundamental energy source for an SEE lies in current- 
carrying magnetic fields. In such a scenario, the free energy available to power the event is the 
excess "non-potential" magnetic energy - the energy above the minimum-energy, potential 
(i.e., current-free) field to which the field can relax. The estimated available nonpotential 
magnetic energies of the active region producing the SEEs are listed in the column labeled 
'Mag' of Tabled] The estimates were made from full- disk line-of-sight (LO S) magnetograms 
obtained from the Miche lson Doppler Image r (MDI; IScherrer et al.lll995l ) on SoHO, using 
the method described by IWelsch et all (120091 ) and outlined below. 



The 2002 July 20 event (Event #5) was located behind the east limb. Although its 
source is therefore uncertain, subsequent active region maps suggest that this event probably 
occurred in AR 10039, the same source region for the 2002 July 23 event (Event #6). 
Consequently, the estimated non-potential magnetic energy for the latter event is a plausible 
estimate for the non-potential magnetic energy in the former, and indeed the use of such 
an estimate is consistent with the procedure used to estimate the non-potential magnetic 
energy content in other near-limb events (e.g., W72° for 2002 February 20 [Event #1] and 
W84° for 2002 April 21 [Event #2]), for which magnetic fields measurements obtained when 
the pertinent active region was near to disk center were used. However, because there is still 
a finite possibility of misidentification of the active region, we have chosen not to use this 
estimated value either in Table [1] or in the pertinent plots of Section |3j 

Numerous efforts have been undertaken to estimate nonpotential magnetic energies in 
active regions near disk center. The methods include: (1) using the magnetic vir ial theorem 
estimates from chromospheric vector magnetograms (IMetcalf et all Il995l 120051 ); (2) semi- 
empirical flux-r ope modeling usin g Ha and EUV images with MDI line-of-sight (LOS) 



magn etograms (IBobra et all 120081 ); and (3) MHD modeling (IMetcalf et al.lll995t iJiao et al. 



19971) and non-potential field extrapolation based upon photospheric vector magnetograms 
Guo et al-liooa ISchriiver et "alll2008l : iThalmann fc Wiegelmarm1l2008[ IThalmann et al.l2008h 



These methods are labor intensive, and uncertainties in their energy estimates are large. For 
example, error bars on virial free energy estimates can exceed the potential magnetic energy. 
Also, there is considerable sca tter in estimates from studies that employ several methods 
to analyze the same data (e.g.. ISchrijver et al.ll2008f ). A couple of generalizations, however, 



can be made. Free energies determined by virial methods matched or exceeded the potential 
field energy, while free energies estimated using other techniques typically amounted to a few 
tens of percent of the potential field energ y Published values for free energies in analytic 
( ISchrijver et al.l 120061 ) and semi-empirical ( IMetcalf et al.ll2008l ) fields meant to model solar 
fields also hover around a few tens of percent of the potential field energy. 
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We estimated the free (i.e., nonpotential) magnetic energies listed in Tabled] to be 30% 
of the potential magnetic energy determined from MDI full-disk LOS magnetograms. This 
is believed to be a conservative estimate but it has the advantage that it can be readily 
determined for most of the events. Some of the events arose from limb active regions, 
for which simultaneous magnetograms are unavailable. Even if vector magnetograms were 
available, uncertainties in free energies would still be large. With published virial free energy 
estimates ranging to a few times the potential energy, it is possible that the true free energy 
could exceed our estimates by a factor of ~10. 

Apart from two cases where flux was clearly emerging near the time of the event (Events 
#29 and #38), we calculated the potential magnetic energies from magnetograms in which 
each event's source active region was near the disk's central meridian, assuming a rigid 
rotation rate of 13° day -1 . This means the energy estimates were sometimes made a few 
days before or after a given event. Fields were assumed to be radial, so each pixel's line- 
of-sight field strength -Blos was divided by the cosine of the heliocentric angle between the 
pixel and the sub-ob servation point, to generate an estimated radial field, Br. Using a 



Mercator projection (IWelsch et all 120091 ). the corrected pixel values were then interpolated 
onto a 2-D plane. Next, the scalar potential \i where B = — V%, was determined using 
a Green's function method. Finally, the magnetic energy Um was estimated by integrating 
(xBr/Sti) over manually-defined cropping windows that contained each active region. Images 
of the magnetograms used, as well as deprojected data with cropping windows, are online at 



http://solarmuri.ssl.berkeley.edu/~welsch/public/meetings/SADOSC24/ 



In several cases in Table [T] (e.g., Events #22 to 25), the same value of the magnetic 
energy is given for adjacent events up to five days apart from the same active region. This 
is because the magnetic energy was estimated from line-of-sight magnetograms taken when 
the active region was close to disk center. These estimates become increasingly unreliable 
as the active region moves away from disk center. Thus, although it is very likely that the 
active region's magnetic fields evolved substantially over the time between events, there is 
no way to reliably quantify these changes from the available magnetograms. This problem 
will be alleviated with the now regularly avail able vector magnetic field measurements from 



the Helioseismic and Magnetic Imager (HMI; IScherrer et all 120121 ) on the Solar Dynamics 



Observatory (SDO), which can be used to estimate the energy in regions located away from 
disk center. 
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3. Comparisons of Energetic Components 

In this section we present comparisons of the energy contents of the various components 
discussed above, through a series of figures each showing logarithmic scatter plots (the "for- 
est") of the energy content of one component versus the energy content of another, for all 
events ("trees") for which data are available for both selected components. The scatter of 
the points around the logarithmic centroid is due both to the true range of energies of the 
selected events and to the often large uncertainties (up to 2.5 orders of magnitude in some 
cases) in the energy estimates of each component. If the uncertainties are random, then the 
centroid location gives an indication of the average ratio of the energies of the two compo- 
nents being plotted, and the scatter of the points about the centroid provides a measure 
of the overall uncertainty in that ratio. Any "outlier" point indicates an anomalous event, 
which could simply identify an unusually large or small event, or which could reveal intrinsic 
differences in the distribution of energies between the different components or some error in 
the energy estimates for the event in question. 

The component energy comparisons are discussed in the following subsections, with 
associated plots given in Figures [1] to [3j All plots have the same four-order-of-magnitude 
range on each axis, so that the degree of spread in a particular energetic component can be 
readily visualized. In each plot, all events that have measured energies for both components 
are shown. The points are indicated by triangles, except for the "outlier" points lying outside 
the 2a ellipse (see below), which are instead labeled by their event numbers (Tabled]) and 
are located at the center of the respective numbers. Events with only upper and/or lower 
limits are generally not shown. However, we have included pertinent data for the behind- 
the-limb Event #5 and for Event #25, for which there is some ambiguity in the CME kinetic 
energy and we have hence used the geometric mean of the two estimates. The logarithmic 
centroid is shown by a bullseye with its X and Y coordinates, calculated using Equation 
( lAlj) . given in the upper left corner of the plot. The three diagonal dotted lines are lines 
of constant ratio R as defined by Equation (JA2]), with R = Y/X = 100%, 10% and 1%. 
These lines each have tick marks showing the overall "size" of the event A = V XY - see 
Equation (IA3j) . The dashed-line ellipse shows the ±2<r locus; the widths of this ellipse 
perpendicular (Equation (IA4I) ) and parallel (Equation flA5j) ) to the lines of constant ratio 
are measures of the 2a spread in the energy ratio, R, and the event size, A, respectively. 
Points outside this ellipse are considered as "outliers" and will be discussed in Section [5j 

Table [3] lists, for each plot, the energetic components involved, the value of the loga- 
rithmic centroid energies and their ratio, the root-mean-square (RMS) spreads in the values 
of the ratio R and the size A. Also, to quantify possible trends of one parameter vs. the 
other, Table [3] lists the Spearman's rank correlation coefficient p, a quantity that measures 
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the correlation between their rank orders (lowest — > highest) of the variables. The formal 
equations used to determine these different parameters are given in Appendix [XI 

3.1. Radiated Energy in the GOES 1 — 8 A Band vs. Radiated Energy from 

SXR-emitting Plasma 

Figure [T^, shows the scatter plot for the radiated energy in the GOES 1 - 8 A band 
vs. the total energy radiated from the hot SXR-emitting plasma. The points are closely 
rank-correlated (p = 0.96) and also cluster very closely in the perpendicular (R) direction, 
showing that the energy radiated in the GOES 1 - 8 A band is a relatively constant fraction 
(R = 0.05) of the total energy radiated from the SXR-emitting plasma. Indeed, a regression 
analysis of the data in Table [1] shows that the ratio (best estimate ±lcr) of the total energy 
radiated by SXR-emitting thermal plasma in the flare to the observed GOES 1 - 8 A flux is 
15.4 ± 0.8. This strong correlation is not surprising since both plotted energy components 
are calculated from the GOES SXR fluxes; the scatter about the trend line arises from the 
differences in the temperatures of the different events. 

3.2. Thermal Radiated Energy vs. Peak Thermal Energy 

Figure [TJd shows the scatter plot of the total energy radiated from hot SXR-emitting 
plasma vs. the peak thermal energy content of that plasma. The relatively tight correlation 
(p = 0.82) between these two components is expected, since both parameters refer to the 
same SXR-emitting plasma. Event #30 (the M6.4 event on 2005 August 25) is the most 
extreme outlier but it is almost equally weak in both energy components. On average, the 
total energy radiated exceeds the peak thermal energy content by a factor of ~3 (R = 2.7 
in Table [3]), implying continuous re-energization of the SXR-emitting material as the flare 
progresses. 

3.3. Peak Thermal Energy vs. Energy in Flare- Accelerated Nonthermal 

Particles 

Figure [lb shows the scatter plot for the peak thermal energy in the hot SXR-emitting 
plasma vs. the energy in flare-accelerated nonthermal particles (electrons plus ions, when 
available). There is substantially greater spread in the points compared to Figure [lb, but 
nevertheless a reasonable bunching of the points. The maximum spread of the points is 
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Table 2. TIM and FISM Bolometric Energies (10 30 ergs) 



Event 


Date 






TTM 




FISM 


Difference 


No. 




Total 1 


Uncertainty 1 


Revised estimate 


Corrected 2 


Corrected 2 


(TIM - FISM) /TIM 


12 


28-Oct-2003 


600 


39% 


362 


362 


310 


14% 


13 


29-Oct-2003 


240 


86% 


137 


137 


128 


7% 


16 


4-Nov-2003 


260 


65% 


142 


426 


447 


-5% 


31 


7-Sep-2005 


300 


71% 


150 


322 


266 


17% 


37 


6-Dec-2006 






46 3 


59 


82 


-39% 



IWoods et all d2006h 

2 Corrected for limb darkening 
3 Moore et al. (2012, in preparation) 



Table 3. Parameters represented in scatter plots. 



Fig. 


Plotted Components 


Log. Centroid* 


R" 


RMSt 




No. of 


No. 


X-axis 


Y-axis 


X 

10 3t) ergs 


Y 

10 30 ergs 




R 


+ 




Events 


la 


Rad. from Hot Plasma 


GOES 1-8 A 


12 


0.6 


0.05 


0.17 


0.51 


0.96 


38 


lb 


Peak Thermal Energy 


Rad. from Hot Plasma 


3.9 


11 


2.7 


0.22 


0.34 


0.82 


26 


lc 


Electrons+Ions 


Peak Thermal Energy 


34 


3.9 


0.11 


0.43 


0.31 


0.36 


26 


Id 


Electrons+Ions 


Rad. from Hot Plasma 


34 


11 


0.31 


0.43 


0.34 


0.46 


26 


2a 


Electrons 


Ions 


32 


11 


0.34 


0.63 


0.52 


0.45 


14 


2b 


CME KE (SW frame) 


SEP 


110 


4.0 


0.04 


0.49 


0.47 


0.47 


20 


2c 


Ions 


SEP 


23 


6.2 


0.27 


0.38 


0.35 


0.20 


6 


2d 


Magnetic 


Bolometric 


890 


55 


0.06 


0.43 


0.38 


0.56 


37 


3 a 


Magnetic 


CME KE+PE 


1000 


200 


0.19 


0.39 


0.34 


0.68 


23 


3b 


CME KE+PE 


Bolometric 


200 


71 


0.35 


0.43 


0.38 


0.54 


21 


3c 


Bolometric 


Electrons+Ions 


49 


34 


0.71 


0.47 


0.35 


0.37 


26 


3d 


Bolometric 


Rad. from Hot Plasma 


57 


12 


0.21 


0.21 


0.47 


0.92 


38 



* X and Y values of the logarithmic centroid, computed using Equation flAH . 

** R = Y/X, the ratio of Y and X values of the logarithmic centroid computed using Equation 1A2I . 

tRMS (root mean square) values of R = Y/X and A = V XY , computed using Equations (I A4t and j A5H ■ The RMS values 
of R and A, respectively, quantify the scatter perpendicular and parallel to the line of constant energy ratio that passes 
through the logarithmic centroid. 

* Spearman's rank correlation coefficient — a non-parametric measure of statistical dependence between two variables - see 
Equation IIA6I I. 
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(d) Thermal radiated energy from the 
SXR-emitting plasma vs. the energy in 
flare-accelerated nonthermal particles (elec- 
trons plus ions when available). 



Fig. 1. — Scatter plots of different energy components, in units of 10 30 ergs. Each plot includes 
all events for which measurements are available for both components. The points are indicated 
by triangles, except for the "outlier" points lying outside the 2a ellipse (see below), which are 
instead labeled by their event numbers (Table [[]) and are located at the center of the respective 
numbers. The location of the logarithmic centroid, defined by Equation f| Al |) . for all the events in 
the plot is shown by a bullseye with its X and Y coordinates listed in the upper left corner of the 
plot. The three diagonal dashed lines represent the 1%, 10% and 100% ratios between the plotted 
components. Lines of constant logarithmic average event energy are shown by dashes every order of 
magnitude along the lines of constant ratio. The major and minor axes of the ellipse are defined by 
±2 times the RMS deviation of the points respectively parallel (Equation (1A5|) ) and perpendicular 
(Equation (IA4I) ) to the line of constant ratio passing through the centroid (see text and Appendix 
for discussion) . Points outside this ellipse are considered as outliers and are discussed in Section [5j 
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(a) Energy in flare-accelerated ions vs. en- 
ergy in flare-accelerated electrons. 



(b) Energy in SEPs vs. CME kinetic energy 
in the rest frame of the solar wind. 
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(c) Energy in SEPs vs. energy in flare- 
accelerated ions. 



(d) Bolometric radiated energy vs. nonpo- 
tcntial magnetic energy in the active region. 



Fig. 2. — Same as Figured] for different combinations of energy components, as indicated on 
the axis labels. 
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(b) Bolometric radiated energy vs. the 
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(c) Energy in flare-accelerated nonthermal 
particles (electrons and ions) vs. bolometric 
radiated energy. 



(d) Total energy radiated by SXR-emitting 
plasma vs. bolometric radiated energy. 



Fig. 3. — Same as Figure [T]for different combinations of energy components, as indicated on 
the axis labels. 
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less than two orders of magnitude in either parameter. On average, the energy in the flare- 
accelerated non-thermal particles exceeds the peak thermal energy by almost an order of 
magnitude (R = 0.11), indicating that there is easily sufficient power in the particles to 
create the SXR-emitting thermal plasma. This conclusion is reinforced by the fact that the 
energy in flare-accelerated electrons is a lower limit (see Section 12 .4[) and is in agreement 
with earlier comparis ons o f flare-accelerated elec t rons v ersus peak thermal energy - see, e.g., 



de Jaeer et all (119861 ) and ISaint-Hilaire & Benzl (120021 ). 



3.4. Thermal Radiated Energy vs. Energy in Flare- Accelerated Nonthermal 

Particles 

Figure [TH shows the scatter plot for the total energy radiated by the SXR-emitting 
thermal plasma vs. the energy in flare-accelerated nonthermal particles (electrons and ions 
when available. This figure combines information already evident in Figure [Tb and Figure [Tb. 
It shows that the energy in accelerated electrons and ions during a flare is not only sufficient 
to supply the peak energy of the SXR-emitting plasma (Figure [Tb), but it is also high enough 
(by a factor of ~ 3, R = 0.31) to account for the radiation from this plasma throughout the 
event. As discussed in Section (2], it follows that a significant fraction of the energy in flare- 
accelerated nonthermal particles is deposited by thermal conduction int o lower-temperatur e 



plasma and ultimately radiated in optical and EUV wavebands (see lEmslie et al.l 120051 ). 
Again, Event #30 (the M6.4 event on 2005 August 25) is the only "outlier" in this plot, 
reflecting the low values of both the thermal and nonthermal energy components. 



3.5. Flare-Accelerated Ions vs. Electrons 

Figure [2h shows the scatter plot for the energy in flare-accelerated ions, as determined 
from the RHESSI gamma-ray observations (Section 12. 5p . vs. the energy in flare-accelerated 
electrons, as determined from RHESSI hard X-ray observations (Section 12. 4p . 

As pointed out in Section I2.4[ the energy in electrons is critically dependent on the low- 
energy cutoff, E min , that is assumed for the electron spectrum. Since the largest value of E min 
that gives an acceptable fit to the data is used for each spectrum, the total electron energy 
values are lower limits with order-of-magnitude uncertainties. As explained in Section I2.5[ 
the situation for the ion energies is even worse, both because of the spread in the observed 
2.223 MeV fluences and because of the need to extrapolate the ion flux at energies above 
30 MeV, as derived from these 2.223 MeV line fluences, to the ion flux above 1 MeV. Because 
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of these large uncertainties in both the electron and ion energies, there is a much wider scatter 
than in the plots in Figure [TJ However, with some notable exceptions that have almost two 
orders of magnitude more energy in the electrons than in the ions (Events #9 and #15), 
the electron and ion energies are genera lly comparable withi n an order of magnitude. This 
result is in agreement with the claims by lRamaty et al.l (119951 ) and lRamaty fc Mandzhavidze 
( 120001 ) and has significant consequences for particle acceleration models. 



3.6. SEP Energy vs. CME Kinetic Energy in the Solar Wind Rest Frame 



Figure |5b shows the scatter plot for the energy in the accelerated SEP population vs. 
the kinetic energy of the CME in the rest frame of the solar wind. We use the solar wind 
rest frame since a shock can be form ed and SEPs accelerate d only if the CME is traveling at 
least as fast as the solar wind speed ( IMewaldt et al.ll2008bl ). Lacking knowledge of the solar 
wind speed low in the corona for each event, we have simply subtracted 400 km s -1 from 
the measured CME speed in order to estimate the kinetic energy available for accelerating 
particles via shock acceleration. 

Most of the SEP values cluster between 1% an d 10% of the CME kine tic energy. Com- 
paring the nine events that are common to both iMewaldt et al.l (l2008al ) and this study 



(Events #2, 4, 7, 8, 11, 12, 13, 14, and 16), the SEP/CME ratio was 5.8% in lMewaldt et al. 



fl2008al ) and is 4% here (R = 0.04 in Table [3}. Overall, result of several changes in the 



analysis, the SEP energy estimates in this paper are reduced from those in IMewaldt et al. 



( 12 0051 . l2008al ) by an average o f 40% . One of these changes is the adoption of the longi- 
tude correction of lLario et al.l ( 120061 ) , which results in changes of as much as a factor of 
two in the ene rgy estimates of individ ual SEP events. For the nine events in common with 
this paper and IMewaldt et al.l (j2008al ). the differences due to longitude/latitude corrections 
alone ranged from -52% to +51% with a mean difference of -16%; for all 20 events shown 
in Figure [2b the average effect is a ~10% decrease in energy. Another change is the adop- 
tion of new corrections for particles crossing 1 AU multiple times an d for adiabatic energ y 
loss that are based on the energy and species-dependent simulations of IChollet et al.l ( 120101 ) : 
these changes resulted in a furth er decrease averaging ~ 30% in the SEP energy content of 
the nine events in common with IMewaldt et al.l ( l2008af ). Finally, the SEP spectra in this 
study were integrated from 0.03 MeV nucieon -1 to 300 MeV nucieon -1 , rather than from 



0.01 MeV nu cieon -1 to 1000 MeV nucieon -1 as in Emslie et al.l (J2004J) and IMewaldt et al. 



( 120051 . l2008al ). The increase in the low-energy limit to 0.03 MeV nucieon 1 re presents a mor e 
realistic threshold for injection into the shock acceleration process (see, e.g., iLi et al.ll2012f ). 



and results in a typical reduction in the energy content by ~5%. The change in the upper 



limit has a negligible effect. 



Overall, this new analysis confirms that the SEP energy is a small, but not insignificant, 
fraction of the CME kinetic energy in most large events. 



3.7. SEP Energy vs. Energy in Flare- Accelerated Ions 

Figure |2b shows the scatter plot for the total energy in SEPs vs. the energy in flare- 
accelerated ions, as determined from RHESSI gamma-ray observations. There are only a 
limited number of events that can be compared but in those few cases there is comparable 
energy in the ions and SEPs (R = 0.27). At first sight, this appears to conflict with the 
study by Mewaldt (2012, in preparation) that shows the number of >30 MeV SEP protons 
measured in interplanetary space is generall y much higher t han the number of >30 MeV 



protons interacting in the solar atmosphere ( jShih et al.l 120091 ) . However, this difference can 



be accounted for in a number of ways. First, as noted in Section 12.51 the inferred energy 
in flare-accelerated ions is based on a very uncertain spectral extrapolation over more than 
an order of magnitude in proton energy (from 30 MeV down to 1 MeV). To obtain the 
energy in flare-accelerated ions (Section 12. 5p . we have assumed a spectral index of 4, which 
results in a significantly higher total energy than would be obtained if the extrapolation was 
performed with the much lower spectral indices representative of SEP spectra measured in 
situ at 1 AU. Further, for the well-observed 2003 October 28 flare (Event #12), we obtained 
an ion spectral index of 3.4 at energies between ~3 and 50 MeV, so that the energy content 
in ions could be significantly lower than we have used here, especially if the spectrum hardens 
even more between, say, 1 and 20 MeV, a feature that is clearly seen in SEP spectra. Second, 
SEP protons typically carry > 80% of the SEP ion energy, whereas flare-accelerated protons 
carry only about one-third of the ion energy. 



3.8. Bolometric Radiated Energy vs. Magnetic Energy 



Figure |2U shows the scatter plot for the bolometric radiated energy vs. the nonpoten- 
tial (free) magnetic energy in the active region. Here we have a relatively tight bunching 
with little more than 1.5 orders of magnitude range in each parameter. It should be em- 
phasized that the plotted bolometric radiated energies are, with the exception of the five 
events noted in Table [p not directly measured but rather estimates made using the FISM 
model (jChamberlin et al.l 120071 120081 ). and the magnetic energy is only good to an order of 
magnitude. With this proviso, we find that the average bolometric radiated energy is ~6% of 
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the free magnetic energy (R = 0.06), and in all cases the available magnetic energy exceeds 
the bolometric radiated energy by at least half an order of magnitude. This is consistent 
with the well-accepted notion that the reservoir of magnetic energy is sufficient to power the 
main components of the flare. 



3.9. CME Energy vs. Magnetic Energy 

Figure [3^ shows the scatter plot for the CME total energy (potential + kinetic) in the 
rest frame of the Sun vs. the nonpotential energy in the magnetic field. The CME energy 
is, on average, only a small fraction (R = 0.19) of the magnetic energy, similar to the value 
R = 0.06 found for the ratio of bolometric radiated energy to magnetic energy (Section |3.8( 
Figure [2H). While bearing in mind the very approximate values of the latter, it nevertheless 
appears, from the results of this and the previous subsection, that much of the available 
magnetic energy (some two-thirds) is retained in the active region (i.e., the field does not 
return to a fully potential state), even after the flare and the ejection of t he CME. This result 



is consistent with the generous limits established by iMoore et al.l (120121 ) on the possible free 



energy that an active region magnetic field can hold before it erupts. 



3.10. Bolometric Radiated Energy vs. CME Energy 

Figure EH and show, respectively, that the available magnetic energy is about 15 
times the bolometrically radiated energy in the flare (R = 0.06) and about 5 times the CME 
energy (R = 0.19). Given the substantial overlap of events common to both plots (Tables [1] 
and [3]), this indicates that, on average, the energy in the CME is larger than the energy 
radiated by a factor of ~3. Figure [3b confirms this result by showing the scatter plot for 
bolometric radiated energy vs. the CME total energy (kinetic + potential). On average, the 
bolometric energy is indeed about half an order of magnitude less than the CME energy 
(# = 0.35). 



3.11. Flare- Accelerated Particle Energy vs. Bolometric Radiated Energy 

Figure [3fc shows the scatter plot for the total energy in flare-accelerated particles (elec- 
trons plus ions) vs. the bolometric radiated energy. This figure shows that the energy in 
accelerated particles during a flare is comparable to the total bolometric radiated energy 
from the flare (R = 0.71), with the ratio being greater than unity in some events and less 



-29 - 



than unity in others. It must be recalled that while the bolometric radiated energies are ac- 
curate to within a factor of ~2 to 3, the energies in flare-accelerated particles are uncertain 
to at least an order of magnitude. The energies in electrons are most probably lower limits 
and may well underestimate the true energy content by up to an order of magnitude. The 
energies in ions may, however, be overestimates, depending on the power-law index used for 
the spectral extrapolation. We can tentatively conclude, however, that there is sufficient 
energy in the flare-accelerated particles to account for all the energy radiated in the flare. 
However, this conclusion must somehow be verified by more accurate estimates of the energy 
in the flare-accelerated electrons and ions. 



3.12. Energy Radiated by Soft X-Ray Emitting Plasma vs. Bolometric Energy 

Figure |3H shows the scatter plot for the total energy radiated by SXR-emitting plasma 
vs. the bolometric radiated energy. The relatively tight correlation (p = 0.92) is to a large 
extent due to the fact that most of the bolometric radiant energies were computed using 
the FISM model, which uses the radiated energy from SXR-emitting plasma to estimate the 
bolometric energy. Nevertheless, the data show that only about one-fifth of the bolometric 
energy radiated by a flare is radiated by the SXR-emitting plasma (R = 0.21). 



Discussion 



The comparisons of the energetics of the different components shown in the scatter plots 
of Figures [TH3] are summarized in Table [3j 

In attempting to draw any definitive conclusions from these comparisons, we must keep 
in mind the limitations of our analysis. The events in our list cover the period from RHESSFs 
launch in February 2002 through 2006. They include nine of the eleven X5 or greater events 
that occurred in this time period. They also include the largest SEP events and all significant 
RHESSI gamma-ray line events, and hence represent events where a significant number of 
ions were accelerated to high energies, either in the flare, at the CME shock, or at both 
locations. Most of the C MEs that we have i nclud ed have kinetic energies of > 10 32 erg. 



According to Figure 8 of iGopalswamy et al.l (120041 ). only 16 CMEs observed from 1997 - 



2002 had kinetic energy >1.4xl0 ergs. Adding in the five CMEs with kinetic energy 
>1.4xl0 32 ergs that occurred in 2003 (Events #11, 12, 13, 14, and 16 in Tabled]), we 
see that on ly 21 CMEs with k inetic energy > 10 32 erg occurred during 1997-2003. For 
this period, IGopalswamyi (120061 ) reports a total of 4133 CMEs, with average kinetic energy 
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5 x 10 erg. Therefore, if we consider the 1997-2003 period as typical, only 21/4133 ~0.5% 



o 



of all CMEs have kinetic energies as large as those considered in this paper. In summary, the 
events studied represent the largest solar eruptive events that occurred during the period in 
question. 

Despite the relatively large number of events in our list, the range of energies is typically 
only about two orders of magnitude in any of the energy components. The uncertainties 
on the energy estimates of each component are generally at least an order of magnitude, 
except for the few cases where the bolometric radiated energy is measured with an accuracy 
of a factor of ~2. Thus, in most cases, we cannot expect to see any significant trends 
with the size of the events over the limited range of our selected events. In spite of these 
selection effects and measurement limitations, the scatter plots nevertheless reveal several 
useful results regarding large SEEs in general and about specific events in particular. With 
a few notable exceptions, the points in each plot are bunched together within the expected 
order-of-magnitude uncertainties. However, there is substantially more scatter of the points 
in some parameters than others. Only a few events stand out as outliers in certain plots; 
these outliers are discussed in Section 

The general bunching of the data points in each scatter plot is characterized by the 
logarithmic RMS deviations parallel and perpendicular to the line of constant ratio that 
passes through the logarithmic centroid. The X and Y values of the centroids are shown 
on each plot and listed along with the RMS values in Table The fact that the data 
points generally bunch together can be in terprete d as a n extension to SEEs of the "big 



flare syndrome" (BFS), a phrase coined by iKahlerl (119821 ) based on the strong correlations 



between proton fluxes and associated microwave and hard X-ray burst parameters, and a 
concept which has since come to mean that each flare component scales roughly linearly 
with some absolute measure of fla re "size." At that time, before the so-called "solar flare 



myth" was exposed (lGoslingill993l ). it was not clear that SEPs were generally more likely to 
be accelerated at CME shock fronts, but now the same concept can be applied to include 
CMEs with flares and SEPs. From our data, it is clear that, with some caveats, the BFS 
concept can be applied to all energetic SEE phenomena, including the CME energy. 

Notable exceptions to the scatter of the points being consistent with the expected un- 
certainties in each parameter are the following: 

• The plot of the GOES 1 - 8 A integrated energy vs. the total energy radiated from 
the SXR-emitting plasma over all wavelengths (Figure [Qa,) shows that these two pa- 
rameters are well correlated (the Spearman's rank correlation coefficient p given in 
Table [3] is 0.96, the highest for any pair of parameters). This is not surprising since 
the measurements of the two parameters are not independent - they both use the 
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same GOES X-ray data. The scatter of the points thus reflects only the range of flare 
sizes and the different temperatures - the scatter perpendicular to the line of constant 
ratio, only about half an order of magnitude, is the result of different temperatures, 
while the scatter parallel to that line is almost two orders of magnitude, reflecting the 
range of flare intensities. A regression analysis of the data in Table [1] shows that the 
ratio (best estimate ± standard error) of the total energy radiated by SXR-emitting 
thermal plasma in the flare to the observed GOES 1 - 8 A flux is 15.4 ± 0.8, i.e., that 
the energy radiated in the GOES 1 - 8 A band is about one-fifteenth to one-twentieth 
(R = 0.05) of the total energy radiated by the SXR-emitting plasma over all wave- 
lengths. The presence of the outlier points for Events #3 and 8, several RMS values 
away from the line of constant ratio, suggests that these two events are different in 
that their temp erature is lower than the a verage. This is consistent with the conclu- 



sion reached by iFeldman et al.l ( 119961 ) and iGarcial ( 120041 ) that lower temperatures are 



generally associated with smaller X-ray peaks. 

The range of the ion energies - about three orders of magnitude - is larger than for 
all other parameters. This range is especially evident in Figure E^,, which shows the 
scatter plot for the energy in flare-accelerated ions as determined from the RHESSI 
gamma-ray observations (Section 12. 5j) vs. the energy in flare-accelerated electrons as 
determined from RHESSI hard X-ray observations (Section 12.4ft . Since the ion energy 
content above 1 MeV was deduced by applying the same 5 = 4 spectral extrapolation 
to the energy content above 30 MeV in all events, this large range in energy contents 
above 1 MeV mirrors exactly the spread in the energy contents above 30 MeV and hence 
the spread in the observed 2.223 MeV line fluences. However, it must be noted that 
uncertainty in the value of 5 results in a further large uncertainty in the ion energy 
content above 1 MeV for any specific event, since this quantity is derived through 
spectral extrapolation over one-and-a-half orders of magnitude in ion energy (Section 
2.5). This uncertainty would act to increase the scatter in the energy content above 
1 MeV if the value of 6 was positively correlated with the value of the >30 MeV energy 
content; alternatively, it would act to decrease the scatter in energy content >1 MeV if 
the value of S was negatively correlated with the value of t he >30 MeV energy content. 



In this context, it should be noted that IShih et al.l (120091 ) found a strong correlation 



over more than three orders of magnitude, between the energy in flare-accelerated 
ions and the energy in electrons above 300 keV (rather than the ~20 keV lower cutoff 
energy used here). This suggests that steep [shallow] spectra (high [low] values of 5) are 
associated with low [high] values of the energy content at high energies, and therefore 
that use of individual spectral indices to create more accurate energy estimates of the 
ion energy content above 1 MeV might reduce the scatter in the plot. 
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5. Discussion of Outlier Data Points 

While most of the events in Table [1] lie (by definition) within the 2a ellipses in the 
various cross-correlation plots, there are a few notable exceptions. We now discuss these 
"outliers" and the possible reasons for their unusual energetic partitioning. 

• Event #1. This M5.1 event, on 2002 February 20, is one of the two events with a 
relatively low ratio of CME energy (kinetic + potential) to magnetic energy (Figure^). 
Table [T] shows a paucity of data for other energetic components. Data for GOES 1 - 
8 A emission and total radiated energy from the SXR-emitting plasma are available; 
the event does fall outside the 2a ellipse in Figure [T^i, but only by virtue of its overall 
weakness, not the ratio of GOES 1 - 8 A emission to total SXR-emitting energy. We 
believe that this is therefore simply a weak event, in which only a small fraction of the 
available magnetic energy was dissipated. 

• Event #3. This C5.0 event, on 2002 May 22, has a relatively low ( < 1%) ratio of 
GOES 1 - 8 A emission to total SXR-emitting energy (Figure [T^i). An event with 
limited overall information (Table [TJ, it does, however, appear in Figure [2b, where a 
normal ratio of SEP to CME energy is evident. It should be noted that this is the only 
GOES C-class event in Tabled] and we therefore simply categorize this event as a weak 
GOES event, possibly due to the low temperature (< 9 MK) of the soft-X-ray-emitting 
plasma. This event is the third in a sequence of events starting with an Ml. 5 event 
peaking at 21:29 UT on May 21, followed by a C9.7 event at 00:30 on May 22 and the 
C5.0 event in question at 03:34 UT. RHES 'SI saw parts of each of these events and the 
6-12 keV images show that they came from three distinctly different locations with the 
following spatial centroid coordinates (in arcseconds): (-550, 270), (880, -330), (750, 
-350), respectively. The peak temperatures derived from the GOES data for the three 
events are 13, 11, and 8 MK, respectively. Thus, this event was much cooler than the 
other larger events in our list. 

• Event #8. Data for this M4.6 event, on 2002 November 9, is available for all energetic 
components other than flare-accelerated ions ( Table [T]); consequently, the event appears 
in most of the scatterplots in Figures [TH3J In most of the plots, the event is situated 
within the general bunching of points. However, Figure |3b reveals a low ratio of 
bolometric energy to CME energy, while Figure shows a similarly high ratio of 
energy in accelerated particles to bolometric radiated energy. Together, these point 
simply to an event with a relatively low bolometric radiance, as inferred from the 
FISM model (Section 12.2)) . This arises because of the relatively short duration of the 
GOES event (13 minutes from peak 1 - 8 A flux to 50% of peak). As pointed out by 
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Chamberlin et al.l (120121 ) . "the total radiated output of flares depends more on the flare 
duration than the typical GOES X-ray peak magnitude classification." The relatively 
low temperature of this flare (peak value of 13 MK) also explains the relatively low 
ratio of GOES 1 - 8 A integrated flux to thermal radiated energy (Figure dh). This 
event also appears in Figure [2]d as having a marginally low ratio of SEP energy to 
CME kinetic energy, but this is a presumably unrelated phenomenon. 

Event #9. This XI. 4 event, on 2003 May 27, has the lowest ratio of flare-accelerated 
ion energy to flare-accelerated electron energy (Figured). The ratios of peak thermal 
energy and broad-band SXR radiated energy to energy in flare-accelerated nonthermal 
particles are average (Figures [Tb and [TU), an d the ratio of GO ES 1 - 8 A to total 



SXR-emitting energy is nominal (Figure [Th). IShih et al.l (120091 ) show that this flare 
has a ratio of accelerated >20 MeV nucleon -1 ions to accelerated relativistic electrons 
that is comparable to other gamma-ray flares, suggesting that this flare is an outlier 
because it has a relative deficiency in higher-energy particles, both ions and electrons. 
In particular, the stated ion energy content may be a significant underestimate if the 
ion spectrum is steeper than the 5 = 4 power-law assumed in the extrapolation from 
the >30 MeV proton energy value that is obtained from the observed neutron-capture 
line fluence. 

Event #12. This X17 event, on 2003 October 28, has relatively high CME and magnetic 
energies (Figured and TableUJ). It was one of three events from the "Halloween" active 
region of October-November 2003, in which a very high non-potential magnetic energy 
value (4 x 10 33 ergs; see Figure [2U) was inferred. This event is therefore an "outlier" 
simply because it was a very large event; there are no particularly unusual ratios of 
energetic components. 

Event #17. The data set for this XI. 6 event, on 2004 July 15, is quite extensive, with 
all components measured except CME and SEP energies (which are understandably 
absent given the ~E45° location of the flare). The event has a slightly low ratio of 
bolometric radiant energy to magnetic energy (Figure [2U), and slightly high ratios of to 
flare-accelerated particle energy to bolometric radiated energy (Figure [3b ) and energy 
radiated by SXR-emitting plasma to bolometric radiated energy (Figure EH). It has 
the equal- lowest bolometrically radiated energy content, and the lowest thermal energy 
content of SXR-emitting plasma, of any event studied (Table [1]), but interestingly does 
not show as an outlier in any other plots in which it appears, notably the plot of thermal 
energy content vs. thermal radiated energy (Figure [Tb) and the plot of thermal energy 
content vs. accelerated particle energy (Figure [Tb) (although it is barely inside the 
2cj ellipse in both of these plots). As with Event #8, these event properties may be 
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explained by the relative short SXR duration of this flare - only 4 minutes from the 
GOES 1 - 8 A peak to 50% of peak flux. Although the peak plasma temperature derived 
from the GOES fluxes was 22 MK, the temperature stayed above 15 MK for only about 
8 minutes. The short duration and low temperatures lead the FISM modeling process 
(Section 12.21) to assign a correspondingly low estimate of the bolometric radiant energy 



flChamberlin et al.ll2012f ). 



Event #18. Data for this M7.1 event, on 2004 July 25, is rather limited (Tabled). 
The event appears as an outlier in Figure [2U, due to a low ratio of bolometric energy 
to magnetic energy. Table [T] shows that this event had one of the highest inferred non- 
potential magnetic energy contents of all the events studied, while Figure []Ja shows very 
small values of both radiation in the GOES 1 - 8 A waveband and total energy radiated 
by the SXR-emitting plasma (although the ratio of the two is nominal). RHESSI did 
show a flare flag starting at 05:37 UT, but entered night at 05:42 UT, three minutes 
after the GOES start time; therefore no reliable electron or ion energy measurements 
are available. Further, no CME was observed for this event, which was located at 
approximately W30° longitude. As a result, it is not known whether the event was 
simply very weak radiatively, or whether only a relatively small part of the nonpotential 
magnetic energy available was released. 

Event #25. For this X7.1 event, on 2005 January 20, the ratios of nonthermal particle 
energy to bolometric radiated energy (Figure Eb), soft-X-ray-radiated to peak thermal 
energy (Figure [Tb), and soft-X-ray-radiated to nonthermal particles (Figure [TU) are 
nominal. However, the ion to electron energy ratio (~5:1; Figure [2^ and Table [1]) 
is conspicuously high; indeed, this event had one of the highest ion energy contents 
measured. This high ion-to-electron energy ratio may in part be due to the presence of 
high-energy protons (inferred from the 2.223 MeV line) for several minutes after most 
of the electron-ion bremsstrahlung had dissipated. The event also has a very low ratio 
of CME energy to magnetic energy (Figure Ek), but a very high ratio of SEP energy 
to CME energy (Figure [2b). The ratio of bolometric radiated energy to free magnetic 
energy is nominal (Figure [2H), but it is the only event with a ratio of bolometric 
radiated energy to CME energy that is greater than 100% (Figure [3b). Together, 
these results show that the main component that makes this event an "outlier" is the 
low CME energy. Note, however, that the CME kinetic energy in this event is very 
uncertain. 

The reasons for these unusual circu mstances are not comp letely understood. Much has 



been written about this event (e.g.. iGrechnev et al.ll2008l . and references therein) but 



the acceleration of the intense flux of SEPs has still not been fully resolved. It has been 
suggested that in this event the SEPs were accelerated at the flare site rather than in 
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the CME shock (Ei3 l20Q5al fbh . In similar vein, Isimnettl feood . l2007h concluded that, 
"the relativistic protons were not accelerated by the CME-driven shock." The event 
produced a cosmic-ray ground-level enha ncement that is amo ng the largest recorded in 
the history of cosmic ray measurements (IMishev et al.ll201ll ). How ever, based on the 
partic le spectrum measured over a wide rigidity range (1-20 G V) . iMorgan fc Lopate 
(120081 ) state that, "The 2005 January 20 GLE was an unusual event in its intensity and 
brevity, placing it on the outer edges of parameter space for shock acceleration to GeV 
energies, but still not req uiring a different process, i.e., direct solar-flare acceleration." 
Similarly, iReamesI ( 120091 ) finds that SEPs in this event, as in other ground-level events 
(GLEs), are produced after the shock onset and don' t requi re a separate non-shock 
injection. On the other hand, iMoraal &: McCrackenl (120111 ) conclude that neutron 
monitor observations of this event indicate two separate pulses of high-energy particles, 
one accelerated by the flare, and a second accelerated by the CME-driven shock. 

Event #30. This M6.4 event, on 2005 August 25, had a very low peak thermal energy 
content and a relatively low ratio of the total SXR-emitted energy to the peak thermal 
energy content (Figure [lb). It also has a low ratio of total SXR-emitted energy to 
nonthermal particle energy (Figure [TH). No CME or SEP data are available for this 
East-limb event. For this event, the GOES flux decayed from the M6.4 peak to the C6 
level (10% of the peak flux) in only ~14 minutes, so that the low values of the radiated 
energy in SXR-emitting plasma are due to this simply being a relatively short-lived 
event. 



6. Conclusions 

Despite the rather large uncertainties in the individual measurements used in this anal- 
ysis, the relatively large number (38) of events nevertheless allows us to reach some general 
conclusions about the "typical" ratios of various energetic components in large SEEs. We 
have found the following general statements to hold: 

• Figure [lb shows that the total energy radiated by the SXR-emitting plasma over 
the course of the event exceeds, by about half an order of magnitude (R = 2.8 in 
Table [3]), the peak energy content of the therma l plasm a that produces this radiation. 



This reinforces the conclusions of iMoore et al.l ( 119801 ) that some form of energy is 



continuously supplied to this hot plasma throughout the event; 

Figures [TH and[3fc show that the energy content in flare-accelerated particles (electrons 
and ions) is sufficient to create not only the total energy radiated by the SXR-emitting 
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Fig. 4. — Bar chart showing the (logarithmic) average energies of the different components for the 
six events for which values were obtained for all components - Events #13, 14, 20, 23, 25, and 38. 
The short thin bars show the ±la logarithmic scatter of the energies of the six events. 
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plasma, but also the total bolometric radiated energy of the event; 



Figure shows that the energy contents of flare-accelerated ions and electrons are 
comparable at the or der-o f-magnitude level. This re s ult su pports the earlier claims of 
Ramaty et al.l (119951 ) and iRamaty &: Mandzhavidzd ( 120001 ) and has significant conse- 
quences for acceleration models; 

Figure |2b shows that the SEP energy is typically a few percent (R = 0.04 in Table [3]) 
of the CME kinetic energy in the solar wind rest frame, a result with implications for 
shock-acceleration models of interplanetary particles; 

The combination of Figures |2H, Ek>, and [3fc shows that the available magnetic energy 
is indeed sufficient to power the thermal plasma, flare-accelerated particles, and the 
CME. Although some "double-counting" may be involved in summing these energy 
components (e.g., both the flare-accelerated pa rticles and the CME may transfer energy 
to the ambient plasma; see lEmslie et al.ll2005l ). this result nevertheless conforms to the 
widely-held view that the source of the energy released in SEEs lies in stressed magnetic 
fields. 



Figure H] shows the logarithmic average (i.e., the geometric mean) of the energies with 
±1cj logarithmic scatter of the various energy components for the six events (events #13, 
14, 20, 23, 25, 38) for which all energetic components were measured. (Events #6 and #12 
were not included, since for these events some of the components were determined only as 
upper or lower limits - see Table [TJ) This figure, coupled with the overall ratios summarized 
in Table [3j succinctly demonstrates how, in very approximate terms, the available magnetic 
energy gets distributed in a "typical" flare in our sample: 



1. Of the ~10 33 ergs of available non-potential magnetic energy, approximately 30% is 
released in the SEE, with the remainder staying in the active region as stored magnetic 
energy. Of the ~30% that is released, some 80% (~25% of the available energy) is 
released in the CME (mostly as kinetic energy) and approximately 20% (~5% of the 
available energy) is released as flare-accelerated particles, roughly evenly distributed 
between electrons and ions. 

2. All of the energy in the flare-accelerated particles appears to ultimate ly emerge as 



radia tion across a wide range of wavelengths, from optical to soft X-rays (jEmslie et al. 



20051 ). However, only about one-third of the energy in flare-accelerated particles (~2% 
of the available stored energy) is ultimately radiated from high-temperature soft X-ray- 
emitting plasma. The maximum amount of energy stored as enhanced thermal energy 
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in the soft X-ray-emitting plasma is ~1% of that released, and the amount of energy 
radiated in the diagnostic GOES 1-8 A waveband is only about 5% of the total energy 
radiated by the SXR-emitting plasma, or ~0.1% of the available magnetic energy. 

3. Because of the need for a CME to "overtake" the solar wind and form a shock front 
where SEPs can be accelerated, only about two-thirds of the kinetic energy carried by 
the CME (~15% of the available nonpotential magnetic energy) is available for SEP 
acceleration. The SEP production process is in turn ~4% efficient, so that only about 
half a percent of the released magnetic energy ultimately appears in the form of SEPs. 

Although for completeness we have listed Spearman's rank correlation coefficient (p) 
values in Table [21 little significance can be attached to these values other than for the 
obviously tight correlations between parameters that are essentially derived from the same 
data (e.g., the bolometric emission, the energy radiated from the SXR-emitting plasma, 
and the energy radiated in the GOES 1 - 8 A band, all of which are dependent on GOES 
SXR flux measurements). Any correlations amongst independent components are masked 
by the large uncertainties in the individual measurements used in the various scatter plots. 
Progress in this direction will require sampling of events over a much larger range in flare 
size to determine if the distribution of energies amongst the different components found 
here for large events is preserved for smaller events, as would be expected for a u big SEE 
syndrome." Such a project is the next step towards a more comprehensive understanding of 
energy release in SEEs. 

Using data from the new, more sophisticated, instruments that are now available will 
allow more accurate energy estimates to be made of some of the components. For example, 
the detailed differential emis sion measure analysi s now possible using data from the EUV 
Imaging Spectrometer ( EIS; ICulhane et al.l 120071 ) on Hinode, and from the Atmospheric 
Imaging Asse mbly (AIA; Lemen et al. 2012h and the Extreme Ultraviolet Variability Exper- 
iment (EVE; IWoods et al.ll2012l ) on the Solar Dynamics Observatory (SDO), should lead to 
better estimates of the energies in the thermal plasma. Ve ctor magnetograms from the He- 
lioseismic and Magnetic Imager (HMI; IScherrer et al.ll2012l ) on SDO allow for more accurate 
estimates of the energy in the nonpotential magnetic field. 

Other energy components not considered in our analysis may be found to contain sig- 
nificant total energy and should be included in any future compilation of global energetics. 
These include the turbulent mass motion s revealed by the br oadening of atomic lines seen 
with EUV and X-ray spectro meters (e.g., iPhillips et al.ll2008l ). and the cumulative heating 



of CME plasma reported by iMurphy et al.l ( 1201 ll ) to be comparable to (or even greater 
than) the CME kinetic energy. An other aspect n ot dis cusse d here is t he qu estion of a sec- 
ond flare phase that, according to IWoods et al.l ( 1201 ll ) and ISu et al.l (j201l[ ). can release a 
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similar amount of energy as in the initial phase. Nevertheless, we believe that the order- 
of-magnitude comparisons of energetic components presented herein represents a significant 
advance in our understanding of the nature of energy release in SEEs. 
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A. Appendix 

Here, for definiteness, we provide the equations used to determine the values listed in 
Table [31 the (logarithmic) centroid energies X centroid and Y centroid , their ratio R, the RMS 
(lcr ) values both perpendicular and parallel to the lines of constant ratio, and the Spearman's 
rank correlation coefficient p. 

The coordinates of the logarithmic centroid in each plot are given by 



1°§10 -^centroid 
1°§10 Ycentroid = 

where N is the number of events for which 
included in the scatter plot in question. 

The lines of constant ratios (R = 100%, 
satisfy the following relation: 



= -J2 lo sio^ ; 

i=l 
1 - 

= ^E lo Sio^ > ( A1 ) 

i=l 

there are viable measures of both components 
L0%, and 1%) between the X and Y components 



R = Y/X 
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or \og 10 R = hg w (Y/X) 

= !ogio Y - lo Sio x ■ 



(A2) 



Lines of constant logarithmic average event energy (\og 10 A), shown in the plots as dashes 
every order of magnitude along the lines of constant ratio, are defined as follows: 



log 10 A = (log 10 X + log 10 F)/2 
= log 10 VXY . 



(A3) 



The RMS deviations of the points perpendicular to and parallel to the line of constant ratio 
passing through the centroid are defined by 



RMS i 



( 1Ogl Ri ~~ 1O Sl0 R centroid) 2 

(l/AO ^Yl ( l °gw( Y i/ X i) ~ \og 10 (Y centroid / x centroid )y 



(A4) 



and 



RMSu 



10 ^centroid) 2 



= (VAOyfe^sioA-iog 

= (l/2N)^2 [log 10 (Xili) - \og w (X centroid Y centroid )} 2 . 



(A5) 



These RMS values were used to draw the ellipse with axes of 2xRMSj_ and 2xRMS|| around 
the logarithmic centroids in each plot. 

The Spearman's rank correlation coefficient p is calculated by first assigning ranks Xj (= 
1, . . . , N) and tji (= 1, . . . , N) to the X and Y values, respectively, of the N points used in 
the plot in question. The ranks are assigned such that X(xi) < X(x i+ i), Y(y,i) < Y(y i+1 ), 
with "ties" assigned the average rank of the tied values. (Note that the rank order does not 
depend on whether we use the X^ values or their logarithms log 10 Xj.) The Spearman's rank 
correlation coefficient is then calculated as the correlation coefficient of the ranks: 



P = 



Y, i=1 (xi - x)(yi - y) 



(A6) 
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where x and y are the means of the ranks Xi and respectively. For a monotonic dependence 
of Y on X, p = 1, even if the variables do not obey a perfect linear correlation. 
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